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Abstract 

Amorphous  solid  electrolytes  in  the  system  xLi20-(l  —  x)(0.6Li2S-0.4P2S5)  (with  x  =  0.1,  0.15,  0.2,  0.25,  0.3,  0.35  and  0.4)  have  been  prepared 
from  a  mixture  of  Li20,  Li2S  and  P2S5  using  a  mechanical  milling  technique  at  room  temperature  and  in  inert  atmosphere.  The  oxysulfide  powder 
mechanically  milled  for  80  h  with  x  =  0.25  shows  an  ionic  conductivity  (cr)  of  2.5  x  10  4  S  cm  1  at  300  K  which  is  the  highest  in  this  series.  The 
glass  transition  temperature  ( Tg )  exhibits  a  minimum  at  x  =  0.25.  Constant  volume  molecular  dynamics  was  done  for  all  x.  Structural  and  transport 
results  from  molecular  dynamics  (MD)  simulation  are  described.  Changes  in  structure  with  x  and  their  influence  on  cr  and  To  are  discussed. 

©  2006  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

Very  few  structural  reports  based  on  mechanically  milled 
oxysulfide  glasses  exist  [1-4].  29Si  MAS-NMR  and  XPS 
studies  for  melt  quenched  glass  system  (100  —  x)(0.6Li2 
S-0.4SiS2)-*LLMOv  (LixMOv  =  Li4Si04,  Li3P04,  Li4Ge04, 
Li3B03,  Li3A103,  Li3Ga03  and  Li3In03)  was  reported  by 
Minami  et  al.  [5].  These  studies  reveal  that  in  the  systems  with 
small  amounts  of  LLMOy,  the  silicon  atoms  coordinated  with 
both  bridging  oxygen  atoms  and  non-bridging  sulfur  atoms. 
The  presence  of  large  amounts  of  these  unique  silicon  struc¬ 
tural  units  brought  about  high  conductivities  and  high  thermal 
stability  of  these  glasses.  On  the  other  hand,  presence  of  non¬ 
bridging  oxygen  atoms  and  S2~  units  created  by  large  amounts 
of  LixMOv,  lowered  these  excellent  properties.  Crystalliza¬ 
tion  kinetics  and  structural  changes  during  crystallization  of 
(100  —  x)(0.6Li2S-0.4SiS2)-xLi3PO4  oxysulfide  glasses  were 
investigated  by  Hayashi  et  al.  [6].  The  composition  dependence 
of  the  activation  energy  showed  opposite  variation  to  the  dif¬ 
ference  Tc-Tg,  where  Tg  and  Tc  are  the  glass  transition  and 
crystallization  temperatures,  respectively.  This  inverse  correla¬ 
tion  was  explained  in  terms  of  the  fragility  of  the  glasses  studied. 
Solid  state  MAS-NMR  measurement  revealed  that  silicon  and 
phosphorus  atoms  coordinating  with  both  sulfur  and  oxygen 
atoms  present  in  oxysulfide  glasses  were  not  detected  after  crys¬ 
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tallization.  In  mechanical  milling,  when  Li20  and  P2O5  were 
separately  used  as  oxygen  source  in  Li2S-P2S5  system  [7],  the 
authors  found  that  P2O5  reacted  more  effectively  compared  to 
Li20  in  the  formation  of  glass.  In  order  to  understand  the  short 
range  order  around  P  and  Li  atoms  and  the  emerging  structural 
changes  due  to  the  addition  of  LEO,  molecular  dynamics  study 
ofxLECMd  —  x)(0.6Li2S-0.4P2S5)  with  x  =  0.1, 0.15, 0.2, 0.25, 
0.3,  0.35  and  0.4  was  taken  up  and  reported  here.  To  analyze  the 
MD  generated  configuration  of  atoms,  graph  theory  techniques 
[8,9]  were  used  to  completely  elucidate  the  connectivities  of  the 
various  units  present  and  their  ratios.  Statistics  of  chains  and 
rings  present  were  also  calculated  using  this  technique. 

2.  Glass  preparation  and  characterization 

The  starting  materials,  Li20,  Li2S  and  P4Sio  (Aldrich,  99%) 
were  placed  in  an  agate  pot  of  volume  50  ml  with  10  agate 
balls  of  10  mm  diameter  and  2  agate  balls  of  8  mm  diameter. 
The  mechanical  milling  was  conducted  using  a  high-energy  ball 
mill  at  room  temperature.  The  rotation  speed  of  base  and  bowl 
was  fixed  at  380  rpm.  The  ball  to  mixture  weight  ratios  were 
in  between  10  and  20.  All  the  processes  were  carried  out  in  a 
glove  box  that  was  filled  with  dry  nitrogen  gas.  The  diffrac- 
togram  taken  at  the  start  of  the  milling  process  for  the  glass 
with  x  =  0.2  showed  peaks  due  to  Li20,  Li2S  and  P4Sio  (Fig.  1). 
At  the  end  of  80  h  of  milling  all  the  peaks  disappeared  leav¬ 
ing  a  halo  pattern.  Fig.  2  shows  the  X-ray  diffraction  of  some 
of  the  glasses  after  80  h  of  milling.  The  density  of  these  pow¬ 
ders  were  measured  using  Archimedes  principle  taking  toluene 
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Fig.  1.  XRD  of  0.2Li20-0.8(0.6Li2S-0.4P2S5)  with  different  periods  of 
milling. 

(density  0.863  gem-3)  as  the  immersing  medium  (Table  1). 
Ionic  conductivity  measurement  was  done  at  room  tempera¬ 
ture  by  Keithley  3330  impedance  analyzer  in  the  frequency 
range  40-100  kHz.  Pellets  of  the  sample  of  size  12  mm  (diam¬ 
eter)  x  2  mm  (thickness)  were  made  using  Ag  powder  on  both 
sides  of  the  sample  and  applying  a  pressure  of  3  tonnes.  DC  con¬ 
ductivity  was  calculated  from  the  plot  of  log/ versus  logtx,  by 
extrapolating  the  frequency  independent  plateau  region  to  zero 
frequency.  The  glass  transition  temperature,  Tg  was  obtained 
using  a  NETZSCH  DSC  200  PC,  at  a  heating  rate  of  1 0  °C  min- 1 
under  constant  flow  of  high  pure  N2  gas.  Measured  values  of 
density,  a  and  Tg  are  given  in  Table  1  and  in  Figs.  3  and  4. 

3.  Molecular  dynamics  (MD)  simulation 

Constant  volume  MD  was  performed  on  cubic  cells  calcu¬ 
lated  with  measured  densities  (Table  1).  Starting  at  a  temperature 


Fig.  2.  XRD  of  xLi20-(l  -x)(0.6Li2S-0.4P2S5)  after  80h  of  milling. 


Fig.  3.  Variation  of  measured  Tg  with  x  for  xLEO-Q  —  x)  (0.6Li2S-0.4P2S5) 
system. 


Table  1 

Physical  parameters  forxLECMl  —  x)(0.6Li2S-0.4P2S5)  system 


Parameter 

Density  (gem  3) 

a  (experimental)  (xlO  4Scm  ]) 

Tg  (experimental)  (K) 

0  (MD)  (x  10-2  S  cm-1) 

Tg  (K)  (MD) 

* 

11 

0 

O 

1.86 

0.73 

382.5 

23.55 

476 

* 

II 

O 

Lf\ 

1.89 

0.44 

384.6 

8.24 

485 

x  =  0.20 

1.92 

0.85 

386 

29.59 

508 

x  =  0.25 

1.86 

2.51 

381 

48.31 

467 

x  =  0.30 

1.82 

2.21 

388 

30.44 

520 

x  =  0.35 

1.72 

0.35 

386.6 

3.82 

495 

x  =  0.40 

1.79 

0.63 

386 

43.82 

502 
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Fig.  4.  Variation  of  measured  a  with  x  for  0.2Li20-0.8(0.6Li2S-0.4P2S5)  sys-  Fig-  5.  PCF  and  RCN  of  P-0  for  0.2Li20-0.8(0.6Li2S-0.4P2S5)  at  300  K. 
tem. 


of  3000  K,  the  system  was  cooled  by  scaling  velocities  to  2000  K, 
1500  K,  1000  K,  750  K,  600  K,  350  K  and  finally  to  300  K  with 
time  step  size  of  10“ 15  s.  At  each  temperature,  the  system  was 
equilibrated  for  20  ps  and  various  parameters  such  as  pair  corre¬ 
lation  function  (PCF)  and  running  coordination  number  (RCN) 
for  various  atomic  pairs  were  averaged  for  another  80  ps.  At 
600  K,  350  K  and  300  K,  the  equilibration  and  averaging  times 
were  increased  to  lOOps  each.  The  interatomic  potential  used 
was  of  the  form 


Z/Z.-e 2 

<t>ij  =  - —  +  Aij  exp 


+ 


8 iOjik  -  n)2 


X 


(.00  -  Ttf 


(0 jik  -  It)2 


exp 


f  rij  T  rik\ 
\  PU  ) 


where  Z\  is  the  atomic  number  of  atom  i  and  Go  is  the  equi¬ 
librium  angle  subtended  by  r/y  and  Ay,  pa0 ,  ky  and  py  are 
constants  of  simulation.  The  third  term  in  the  potential  is  the 
three  body  vessel  term  introduced  to  account  for  the  covalent 
nature  of  interactions.  Z;  for  Li,  P,  O  and  S  are  0.4e,  2e,  — 0.8e  and 
— 0.8e,  respectively.  Initial  values  for  constants  of  simulation  and 
atomic  number  Z;  were  obtained  from  MD  studies  reported  for 
the  binary  systems  xLi2S-(l  —  v)P2S5  [10].  They  were  subjected 
to  further  adjustments  till  equilibrium  was  attained  at  all  sim¬ 
ulated  temperatures  with  chemically  acceptable  bond  lengths. 
Table  2  contains  the  relevant  potential  parameters  used. 


Table  2 


Potential  parameters  used  in  two  body  Buckingham  term  for  xLFO-Q  —  x) 
(0.6Li2S-0.4P2S5)  system 


Parameter 

Aij  (eV) 

Pijo  (A) 

P-P 

4000.91 

0.35 

P-O 

190.0 

0.36 

P-S 

1000.57 

0.288 

S-S 

4200.43 

0.27 

Li-S 

1000.42 

0.26 

Li-O 

200.0 

0.28 

To  calculate  the  mean  square  displacement  [MSD]  for  all 
atoms,  the  averaging  time  was  increased  to  lOOps  at  all  tem¬ 
peratures.  The  diffusion  constant  D  for  Li+  was  obtained  from 
the  slope  of  the  linear  part  of  the  mean  square  displacement 
versus  time  interval  (r)  curve.  This  follows  from  the  Einstein’s 
equation, 

D=\jx  “  r(yt°  +  T^2) 

where  ([r(fo)  —  r(to  +  r)]2 )  is  the  MSD,  r  the  time  interval  and 
(•  •  •)  indicates  average  over  time  origins. 

Tg,  the  glass  transition  temperature  was  obtained  from  the 
change  in  slope  in  the  D  (for  Li+)  versus  temperature  plot.  The 
structure  of  the  simulated  glass  was  analysed  from  connectivity 
matrices  and  graph  theory  techniques  as  reported  by  us  earlier 
[11,12]  and  averaged  over  50  different  equilibrated  configura¬ 
tions  at  300  K. 

4.  Results  and  discussion 

PCF  and  RCN  for  P-0  are  shown  in  Fig.  5  at  300  K.  MSD 
versus  time  interval  (r)  for  Li,  P,  O  and  S  are  shown  in  Fig.  6 
for  v  =  0.4  at  300  K.  Such  plots  for  other  glasses  are  similar,  a 
and  To  obtained  from  MD  are  given  in  Table  1 .  Analysis  of  the 
equilibrium  configuration  from  MD  using  the  adjacency  matrix 
yielded  parameters  in  Table  3.  In  this  table,  r;_y  denotes  the 
average  distance  between  atoms  i  and  j  and  rti-j  denotes  the 
average  number  of  j  atoms  around  an  i  atom.  BO  represents 
oxygens  in  P-O-P.  NBS  represents  sulfurs  that  do  not  take  part 
in  P-S-P  bridging,  but  do  take  part  in  P-S-Li  type  of  linkage 
to  varying  extent.  P-Li-Li-P  denotes  linkages  among  various 
P-S-P  chains  mediated  by  Li  as  P-S-Li-S-Li-S-P. 

The  total  coordination  including  O  and  S  around  P  is  close 
to  5.  Around  Li,  it  varies  between  3.5  and  4.77.  Two  peaks  are 
found  in  P-P  PCF  at  3.67  A  and  4.2  A  for  all  glasses.  These 
peaks  arise  from  P-O-P  and  P-S-P  interactions.  They  are  of 
varying  heights  due  to  changing  fractions  of  these  bridges.  All 
distances  r/_y  show  little  variation  with  x,  with  exception  of  rp_o 
which  has  large  values  for  v  =  0. 1  and  0.15.  This  is  also  observed 
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Fig.  6.  MSD  vs.  t  for  0.4Li20-0.6(0.6Li2S-0.4P2S5)  at  300  K. 

in  a  similar  system  vLi20-(l  —  x)P2S5  [13]  which  occurs  for  a 
low  ratio  of  ?zp_o/ftp-s,  the  reasons  for  which  is  not  clear. 

The  coordination  around  P  is  made  up  of  units  such  as  PS5, 
POS4,  PO2S3  and  PS3  in  varying  fractions  for  the  different 
compositions.  The  reported  unit  POS3  at  the  composition  for 
the  highest  sigma  found  in  the  glass  Li2S-P2S5-P205  [7]  also 
occurs  in  this  system  at  v  =  0.20,  which  is  close  to  the  maximum 
conductivity  region  of  v  =  0.25.  For  all  compositions,  nearly 
equal  fractions  of  OS3  and  LiS4  units  are  found  around  Li  with 
minor  fraction  of  L1S5.  For  v  greater  than  0.25,  units  such  as 
LiOS2,  LiOS3  and  LiOS4  units  are  found  in  smaller  numbers. 
Chain  structure  exists,  made  up  of  P-S-P  chains.  The  average 
chain  length  decreases  from  11.4  to  8.8  at  v  =  0.25  and  then 
increases  to  12.9  at  v  =  0.4.  The  average  chain  length  made  up 
of  P-O-P  bridges  is  slightly  above  1  and  varies  between  1.03 
and  1.38  which  shows  that  P-O-P  chains  do  not  exist. 

Tg  variations  with  v  for  measured  and  simulated  data  are  given 
in  Table  1 .  Their  increase  and  decrease  with  v  follow  the  same 


order  except  at  the  end  between  v  =  0.35  and  0.4.  The  fraction  of 
bridging  oxygens  which  strengthens  the  glass  has  a  minimum  at 
v  =  0.25  and  understandably  the  Tg  also  has  a  minimum  there. 
The  rise  and  fall  of  a  with  v  is  similar  for  the  simulated  and  mea¬ 
sured  data  except  at  v  =  0.4,  where  the  simulated  a  shows  a  sharp 
increase.  Simulated  a  is  three  orders  of  magnitude  higher  than 
those  of  the  measured  glasses.  Higher  a  is  generally  observed  in 
simulated  glasses  due  to  high  Active  temperatures  caused  by  the 
very  high  cooling  rate  in  simulation  (109Ks-1).  Ball  milled 
glasses  have  even  lower  cooling  rate  compared  to  quenched 
glasses  which  gives  rise  to  the  large  difference  in  sigma  observed 
here. 

Inspection  of  Table  3  shows  that  many  structural  changes  are 
taking  place  at  v  =  0.25,  in  addition  to  the  maximum  in  a  and 
the  minimum  in  Tg  both  of  which  occur  for  v  =  0.25.  Structural 
changes  observed  at  v  =  0.25  and  their  causes  are  given  below: 

(1)  np_s  remains  high  up  to  0.25  and  then  decreases  and  nus 
becomes  a  maximum.  Three  factors  contribute  to  this  coor¬ 
dination  behaviour.  Both  Li  and  P  share  the  sulfur s,  P/S  ratio 
remains  at  0.307  and  #(Li  +  P)/#S  ratio  becomes  greater  than 
1  for  v>  0.25,  giving  rise  to  lower  ny>s  and  nn-S  at  higher 
v. 

(2)  Both  Li-O-P  and  Li-S-P  have  maximum  values.  The  frac¬ 
tion  of  NBOs  having  more  than  one  Li  around  them  is 
maximum  (0.33(3))  for  v  =  0.25  resulting  in  the  maximum 
for  Li-O-P  interaction.  The  fraction  of  bridging  oxygens 
interacting  with  Li  is  negligibly  small  for  all  v.  But  for  sul¬ 
fur,  the  opposite  happens.  The  fraction  of  BS  having  more 
than  one  Li  around  it  is  maximum  for  v  =  0.25  which  causes 
the  maximum  value  in  Li-S-P. 

(3)  The  average  P-S-P  chain  length  is  a  minimum  even  though 
it  coincides  with  the  minimum  in  NBS  population.  Short 
P-S-P  chains  with  smaller  NBS  population  is  possible  if 
the  chains  possess  many  branches  which  leads  to  higher  Mw 
for  n  less  than  the  chain  length.  For  x  =  0.25 ,  the  average 
chain  length  is  8.8  and  its  M6  is  on  the  higher  side. 


Table  3 


Parameters  obtained  from  graph  theory  analysis  forxLi20-(l  —  x)(0.6Li2S-0.4P2S5)  system 


Parameter 

X 

11 

0 

0 

X 

11 

0 

U\ 

x  =  0.20 

x  =  0.25 

x  =  0.30 

x  =  0.35 

x  =  0.40 

rP_o  (A) 

1.82 

1.82 

1.62 

1.59 

1.62 

1.62 

1.62 

rp-S  (A) 

2.12 

2.12 

2.12 

2.15 

2.15 

2.15 

2.15 

m-o  (A) 

2.07 

2.07 

2.07 

2.06 

2.07 

2.07 

2.07 

ru-s  (A) 

2.52 

2.52 

2.52 

2.52 

2.52 

2.52 

2.52 

ftP-0 

0.214(6) 

0.286(9) 

0.419(7) 

0.493(9) 

0.68(1) 

0.80(1) 

0.93(1) 

ftp-s 

4.65(1) 

4.54(2) 

4.53(2) 

4.66(2) 

4.27(1) 

4.04(2) 

4.03(3) 

ftLi-0 

0.050(4) 

0.109(9) 

0.156(7) 

0.30(1) 

0.26(1) 

0.34(1) 

0.44(1) 

ftLi-S 

3.60(3) 

3.50(3) 

3.63(3) 

4.47(2) 

3.38(3) 

3.18(3) 

3.19(2) 

BO 

0.54(4) 

0.38(2) 

0.36(2) 

0.24(1) 

0.31(2) 

0.25(1) 

0.17(1) 

NBS 

0.449(5) 

0.474(9) 

0.489(5) 

0.451(6) 

0.569(7) 

0.614(1) 

0.609(1) 

s2- 

0.031 

0.046 

0.046 

0.052 

0.048 

0.062 

0.071 

#(Li-0-P)/#0 

1.6(1) 

1.9(1) 

2.2(1) 

3.6(1) 

2.7(1) 

3.0(1) 

3.4(1) 

#(Li-S-P)/#S 

10.6(1) 

10.9(1) 

12.3(1) 

18.7(1) 

13.0(1) 

12.7(2) 

14.4(1) 

Average  (P-S-P)  chain  length 

11.4(3) 

10.3(4) 

10.8(5) 

8.8(3) 

12.8(7) 

12(1) 

12.9(2) 

Average  (P-O-P)  chain  length 

1.17(3) 

1.037(3) 

1.23(2) 

1.14(4) 

1.38(1) 

1.359(9) 

1.49(6) 

S/O 

23.4 

14.73 

10.4 

7.86 

6.06 

4.82 

3.9 

rt-j  values  are  the  peak  values  obtained  from  MD. 
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(4)  Tg  registers  a  minimum  where  the  BO  population  shows  a 
dip  though  it  does  not  have  a  minimum  there.  To  explain  this 
behaviour  of  Tg  at  all  x,  it  is  necessary  to  take  the  oxygen 
content  (OC)  of  the  glass  into  account  as  well.  If  G  denotes 
the  product  of  BO  and  OC,  then  the  behaviour  of  G  with  v 
is  similar  to  that  of  Tg. 

(5)  a  is  maximum.  For  the  compositions  we  have  studied  here, 
the  BO  population  is  high  for  v  less  than  0.25  and  NBS 
population  is  high  for  v  greater  than  0.25.  The  S2~  fraction 
increases  mildly  with  v.  The  first  two  parameters  have  low 
values  for  v  =  0.25  and  yet  a  attains  maximum  value  there. 
However,  for  this  composition,  the  average  number  of  Li+ 
found  near  BO  and  the  corresponding  number  for  BS  have 
maximum  values  of  0.11  and  0.82  among  all  other  com¬ 
positions,  which  creates  a  very  conducive  environment  for 
higher  Li  mobility.  At  v  =  0.4,  the  NBS  population  is  on  the 
higher  side,  though  the  other  factors  are  not  so  favourable. 
Hence  its  high  a  can  only  be  attributed  to  the  large  NBS 
population  found  there.  It  is  pertinent  to  add  here  that  the 
measured  a  does  not  show  such  a  high  value  for  v  =  0.4. 

It  is  seen  that  for  the  ratio  of  sulfur  to  oxygen  content  (S/O) 
between  6  and  7.8,  experimental  a  has  the  highest  range  of 
values  between  2.21  and  2.51  x  10~4Scm_1.  High  values  in 
this  range  are  also  obtained  for  o  from  MD  if  we  neglect  the 
high  a  at  x  =  0A  which  is  not  seen  in  experimental  o.  For 
oxy sulfide  glass  systems  (100  —  v)(0.6Li2S-0.4SiS2)-vLi*MOy 
(Li^MOy  =  Li4Si04,  Li3P04  and  Li4Ge04),  o  is  highest  for  S/O 
equal  to  6.65  and  for  LiAMOv  =  Li3B03,  Li3A103,  Li3Ga03 
and  Li3ln03,  the  highest  a  is  obtained  for  S/O  equal  to  8.8 
[5].  These  values  of  S/O  fall  in  the  range  for  high  a  obtained 
for  our  system.  It  is  surprising  that  for  such  different  systems 
with  varying  coordinations  and  basic  units  in  the  glass  matrix, 
a  seems  to  depend  crucially  on  S/O  ratio.  In  Li2S-P2S5-P205 


glass  prepared  by  mechanical  milling,  the  highest  a  is  obtained 
in  80Li2S-19P2S5-lP2O5,  where  the  S/O  ratio  equals  35.  Higher 
P2O5  content  reduces  a  in  this  system  [7].  This  result  shows  that 
the  source  of  oxygen  matters  in  addition  to  its  fractional  content 
in  the  glass. 

5.  Conclusions 

Both  P  and  Li  coordinate  to  O  as  well  as  S  atoms  in  various 
proportions  in  this  system.  The  BO,  NBS  and  S2_  popula¬ 
tions  control  the  conductivity.  Many  significant  changes  occur 
at  v  =  0.25  where  the  highest  conductivity  and  the  lowest  Tg  are 
observed. 
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